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ABSTRACT 

We show that a cosmic standard ruler can be constructed from the joint measurement of the time delay, At, between gravitationally 
lensed quasar images and the velocity dispersion, tr, of the lensing galaxy. This is specifically shown, for a singular isothermal sphere 
lens, D L K At/it 2 , where D l is the angular diameter distance to the lens. Using MCMC simulations we illustrate the constraints set 
in the £2 m -fi A plane from future observations. 
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1. Introduction 

In the cosmological ACDM model, it is currently estimated 
that 96% of the total energy density of the Universe is 
in the form of dark matter (25%) and dark energy (74%). 
These proportions have been inferred from a number of com- 
pletely independent measurements, e.g. cosmic microwav e 
background (CMB) dSpergel et al.ll2003UKomatsu et al.ll2009l) , 
baryon acoustic oscillations ([E isenstei n et al.ll2005l) . supernovae 
( Riess et al.lll998b | Perlmutter et al.|["l999l) . l arge scale structure 
dPeacocket al.ll200ll) . clusters of galaxies dWhite et all 1 1993b 
and weak gravitational lensing dWeinberg & Kamionkowskil 
120031) . Unfortunately, any astrophysical approach suffers from 
potential systematic uncertainties (e.g. supernovae - not stan- 
dard but standardized candles with possible redshift evolution; 
CMB - various parameter degeneracies and interference with 
foregrounds; weak lensing - PSF influencing the galaxy shape 
measurement). It is therefore important to explore complemen- 
tary methods for measuring these quantities. Gravitationally 
lensed quasars QSOs offer such an attractive alternative. 

QSOs that are positioned in a way that a lens, i.e., a mas- 
sive foreground object such as a galaxy or a group of galaxies, 
intersects the line of sight, can be seen as magnified, multiple im- 
ages. Gravitationally lensed QSOs have already been used to set 
constraints on cosmological parameters. Notably Ikef sdall d 1 9641) 
showed that the Hubble constant Ho can be measured from a 
multiple imaged QSO if the time delay between the lensed im- 
ages and the mass distribution of the lens are known; attempts 
to constrain the cosmological const ant A have been bas ed on 
gravitationally lensed QSO statistics dFukugita et al ] |1990l) . The 
importance of strong gravitational lensing in future constraints 
on the evolution of the dark ene rgy equation of state parameter 
w(z), has also been emphasized dLindej2004l) . 

In recent years, joint studies of stellar dynamics and gravita- 
tional lensing have proven very fruitful, e.g. the Lenses Structure 
& Dynamics (LSD) Survey and the Sloan Lens ACS Survey 
(SLACS) were used to constrain the density profi les of galaxies 
dTreu & Koopmans 1 120041: iKoopmans et al.ll2006l). Methods for 
using either velocity dispersions dGrillo et al. I l2008h or time de- 



lays dDobke et al.ll2009T) of lensed systems have been proposed 
as estimators of Q m and Qa- 

We note that a standard cosmic ruler can be constructed 
from the joint measurement of the time delay between QSO im- 
ages and the velocity dispersion of the lensing galaxy, indepen- 
dent of the redshift of the QSO. We explore its use to constrain 
Q m and Qa in view of the large samples of lenses to be found 
in forthcoming experiments, such as from the Large Synoptic 
Survey Telescope (LSST), Square Kilometre Array (SKA), Joint 
Dark Energy Mission {JDEM), Euclid, and the Observatory 
for Multi-Epoch Gravitational Lens Astrophysics (OMEGA) 
dDobke et al'll2009Hlyezic et al.ll2008fc ICariJli & Rawlingsl l2004 ; 
iMarshall et al.ll2005t iMoustakas et alj|2008l) . 

2. Lenses as standard rulers 

The time delay At is the combined effect of the difference in 
length of the optical path between two images and the gravita- 
tional time dilation of two light rays passing through different 
parts of the lens potential well, 



At 



1 + Zl D os D ol 1 1 



£>LS 



(i) 



Here 6 and fS are the positions of the images and the source re- 
spectively, zl is the lens redshift, and *¥ is the effective gravita- 
tional potential of the lens. £>ol, ^os, ^ls are the angular diam- 
eter distances between observer and lens, observer and source, 
and lens and source, respectively. If the lens geometry 6-/3 and 
the lens potential *P are known, the time delay measures the ra- 
tio DosDol/Dls, also known as the effective lensing distance 
r(z, £l m , Qa), which depends on the cosmological parameters. 

The observed velocity dispersion of a galaxy cr is the re- 
sult of the superposition of many individual stellar spectra, each 
of which has been Doppler shifted because of the random stel- 
lar motions within the galaxy. Therefore, it can be determined 
by analyzing the integrated spectrum of the galaxy, which has 
broadened absorption lines due to the motion of the stars. The 
velocity dispersion is related to the mass through the virial the- 
orem: a 2 oc M lT /R, where M lT is the mass enclosed inside the 
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radius R. The mass is measured by the Einstein angle 6e of the 



lensing system Mg E — 



AC D, 



Ol, where R = Dql^e- Thus, 



~_2 ~. Dps a 

Since the time delay is proportional to Ads^ol/^ls and the 
velocity dispersion is proportional to Dqs/D^s, the ratio At/ a 2 
is dependent only on the lens distance and therefore acts as a cos- 
mic ruler: At /o 2 oc DolH In Fig- 1 we illustrate the dependency 
on the lens redshift of the measurables, time delay At, velocity 
dispersion squared <x 2 , and the ratio between them, At/o 2 . To 
show the sensitivity of the three functions to Qa, we plot them 
for four cases of a fiat (Q m + Qa = 1) Universe with Qa = 
0.2, 0.5, 0.7 and 0.9, relative to an Einstein-de Sitter Universe 

(o m = i,a A = o). 

We can see that At/ct 2 is more sensitive to the cosmologi- 
cal parameters than At or a 2 separately. We also note that the 
higher the lens redshift, the more pronounced is the dependency 
on cosmology, hence it is important for this method to study high 
redshift lenses. Finally, it has the advantage of being independent 
of the source redshift. 



Both velocity dispersion and time delay depend on the grav- 
itational potential of the lens galaxy. For the simple case of a 
singular isothermal sphere (SIS) both parameters can be easily 
expressed analytically. While the SIS model is convenient for its 
simplicity, it is al so a su rprisingly useful model for lens galaxies 
dKoop mans et aUl2006t lOguril l2007t ISchechter & Wambsganssl 
l2004t iGui maraes & Sodre 2 0091: iKoopmans et al.ll2009fc 

The time delay is 



1 + z L DolDqs //l2 

and the velocity dispersion is 
c 2 D os 



°"SIS = 



4tt D ls ' 



(2) 



(3) 



For the SIS model the Einstein angle is half the distance between 
the lensed images, 6e = (0i + 92)12 (we define At and 6*1,2 to be 
positive, with 62 > 6\). Hence 



£>ol«?2 - fli) = 



At 



SIS 



(4) 



3. Monte Carlo Markov Chain simulations 

To explore the constraints set on cosmological parameters from 
the joint measurements of image positions, lens redshift, time 
delay and velocity dispersion, we per form simulations on SIS 
lenses using the Metropolis algorithm (Saha & Will iamslll994l) . 
To simulate the uncertainties of the measurables we use x 1 given 
by: 



X = 



o- 2 fl E (l+--L) 8* 



1 l 2 
- - r(z, Q m , Q A )J 



/cf\ 2 / At SBe \ , ( At Sa 1 \ 2 I SAt Y 



(5) 



1 A standard ruler can be used to measure angular diameter distances. 
Standard candles, on the other hand, measure luminosity distances, 
which can be obtained by multiplying angular diameter distances by 
a factor (1 + z) 2 . 
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Fig. 1. Dependence of the three quantities, At (top panel), a 2 
(middle panel), At/o 2 (bottom panel) on the lens redshift. The 
source redshift was fixed to z = 3 . A flat Universe (Q. m + Qa = 1 ) 
was assumed with 4 different values £1a: 0.2, 0.5, 0.7 and 0.9. 
Each curve is plotted relative to a flat Einstein-de Sitter Universe 
(Q m = 1, Q A = 0). 



where, for simplicity we have assumed that for each simulated 
lens, the lens images are aligned such that 6\ — and 62 = 26^. 
Similar expressions are computed for At and a 2 . In all the sim- 
ulations we let Q m , Qa vary between to 1 while either fixing 
or marginalizing over Hq. 



3. 1 . Comparison of methods 

We first illustrate the constraints set in the Q m -QA plane for 
125 simulated galaxies. Their redshifts are equally distributed 
between 0.1 and 1.1, the source redshifts between 1.5 and 3.5 
and the velocity dispersions in the range 100 to 300 km s~'. For 
each lens we calculate the Einstein angle from Eq. 3 assuming 
Q m = 0.3, Q A = 0.7 and H = 70 km s" 1 Mpc -1 . We con- 
strain the sample to easily detectable systems, thus we include in 
the si mulations only lenses with Qe larger than 0.5" dGrillo et al.l 
l2008h . We assume simulated 5% uncertainties in the Einstein 
angle and run 10000 minimizations of x 2 for three different ob- 
servables: At, o 2 , and At/o 2 (see Fig. 2). 

The constraints on the parameter space from the time 
delay (Fig. 2, top pan el) have a non-linear, curved shape 
(ICoe & Moustakasl2009l) . elongated roughly in the Q m + Qa = 1 
direction. The velocity dispersion (Fig. 2, middle panel) gives a 
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Fig. 2. MCMC simulations for three different cases: At (top, 
panel), o~ 2 (middle panel), and At /a 2 (bottom panel). 125 galax- 
ies were generated with redshifts equally distributed between 0. 1 
and 1.1, source redshifts between 1.5 and 3.5, and velocity dis- 
persions varying from 100 to 300 km s _1 . For each lens we calcu- 
late the Einstein angle for a SIS in a cosmology with Q m = 0.3, 
Qa = 0.7. Hi) was fixed at 70 km s _1 Mpc~'. Yellow regions are 
lcr, green are 2cr confidence levels. 



good constraint on Qa, but a weak one on Q m . Joint observations 
of time delay and velocity dispersion (Fig. 2, bottom panel) give 
constraints approximately perpendicular to the flat Universe line 
similarly to Type la supernovae. As expected, the proposed cos- 
mic ruler (Eq. 4) gives tighter constraints on both cosmological 
parameters than the other two methods. 



3.2. One high-redshift lens 

To study the constraints on cosmological parameters from a sin- 
gle lens we run 10000 minimizations on one syst em with param- 
eters similar to an existing lens, MG 2016+1 12 dLawrence et alj 
11984 . We set the lens redshift to 1.0, the source redshift to 3.27 
and the Einstein angle to 1.7". The velocity dispersion is calcu- 
lated from Eq. 3 assuming Q m = 0.3, Qa = 0.7. We perform 
simulations for two uncertainty scales (5% and 10%). Because 
the results do not depend on the parameter from which the un- 
certainty comes, the simulated error can be understood as the 
uncertainty in either the velocity dispersion squared, the time 
delay, the Einstein angle or a combination of these. We present 
the results from these simulations in Fig. 3 (top row). The prob- 
ability contours form a wide stripe going across the parameter 
space in a direction almost perpendicular to the Q m + Qa = 1 
line. To show the importance of Hq in the simulations we present 
two cases, first for marginalized Ho — 70 + 5 km s _1 Mpc~' (left 
column) and second for fixed Hq = 70 km s _1 Mpc~' (right col- 
umn). The shift in the marginalized case between the 5% and 
10% confidence contours is due to the fact that high values of 
Ho change the distances calculated from the angular diameter 
distance equation less than the low ones, thus, to compensate, 
the region with simulated 10% uncertainties is shifted towards 
lower values of Q m , Qa compared to the 5% region. 



3.3. Many lenses 

We also simulate 20 and 400 lensing systems with lens redshifts 
equally distributed between 0.5 and 1.5 (Fig. 3). In this case 
the probability of detecting the correct cosmological parameters 
converges faster around the assumed 'real' values Q m = 0.3, 
Qa = 0.7, creating, as a constant probability contour, an el- 
lipse. The simulations show that by observing 20-400 lenses 
with small measurement errors, cosmological parameters can be 
well constrained. For a 5% error in the marginalized case we get 
for 20 lenses: Q A = 0.65!° ^, Q m = 0.30!° ^ and for 400 lenses: 



Qa 



Q-7Q+0.05 Q 

U - /U -0.06' " n 



-0.21 

0.30 + o°? 



(1 o~ confidence interval). 



4. Discussion 



In our simulations we have assumed that the lenses have SIS 
mass distributions. The SIS profile seems to be a rather good 
choice, because several studies based on dynamics of stars, 
globular clusters, X-ray halos, etc. have shown that ellipti- 
cal galaxies have approximately flat circular velocity curves 
dKoopmans et alj|2006t lOguril 120071: ISchechter & Wambsganssl 



l2004tlGuimaraes & Sodrel2009l:lGerhard et al.l200ll) . The struc 
ture of these systems may be considered as approximately ho- 
mologous, with the total density distribution (luminous+dark) 
close to that of a singular isothermal sphere. We stress, however, 
that the proposed cosmic ruler does not rely on the assumption 
of a SIS. For any potential characteristic of the lens galaxies it 
is sufficient to assume that the potential (on average) does not 
change with redshift (or that any evolution can be quantified and 
corrected for). In other words, At/ct 2 oc Dol for all lens models. 
While we c annot rule out some redshift evolution, it is expected 
to be weak dHolden et al.ll2.Q09b , given that lens galaxies are typ- 
ically massive and hence fairly relaxed, giving rise to similar 
structures as their lower-redshift counterparts. 

Nevertheless, like for any cosmic ruler, there is a range of 
possible systematic uncertainties which must be addressed be- 
fore it can be used for precision cosmology. Several factors af- 
fect the lensing configuration, the velocity dispersion, and the 
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Fig. 3. MCMC simulations for Ar/cr 2 . Top row: One lens with 
parameters similar to MG 2016+112 (z L = 1.0, z s = 3.27, 
0e = 1.7"). Middle row: 10 lensing systems with lens redshifts 
equally distributed between 0.5 and 1.5. Bottom row: 400 lens- 
ing systems with the same redshifts distribution as for the mid- 
dle row. 10000 MCMC minimizations were performed. All three 
simulations were performed for two uncertainty scales: 5% (yel- 
low) and 10% (green). Left column: The plots marginalized over 
Hq. Right column: Hq is fixed to 70 km s Mpc~' . 



time delays t o various degrees: v elocity anisotropy, total mass- 
profile shape dSchwab et al.l2009h and the detailed density struc- 
tures (mas s profiles of d ark and luminous matter, ellipticities) of 
the lenses (Tonrvfl983l) . mass along the line of sight to the QSO 
dLieul2008l) a nd in the envir onment of the lenses, such as g roups 
and c l usters dMetcalfll2005|) (the mass-sheet degeneracy (jSahal 



I2000t iFalco et all 11985b IWilliams & Sahall200Ct : iQguril 120071)) 



and substructure dPalal & Kochanekll2002f IMaccio & Miranda 



120061: IMaccio 112006b IXull2009ir These factors are of the order of 
our assumed measurement uncertainties. 

Fortunately, the velocity anisotropy ap p ears to be small 



in lens galaxies dKoopmans et al.1 120091: van de Ven et al 



2003) and in relaxed ga laxies in general dHansen & Moore 



2006? lAn & Evansl 120061) . The density structures of ellipti 



tant but usually limited dHolder & Schechterll2003l) and can be 
spotted from anomalous flux ratios or lensing configurations. 

To date, there are around 200 known strong gravitational lens 
systems, but only 20 of them have measured time delays, with 
typical errors of 5-10%. Similar errors are obtained for mea- 
sured velocity dispersions. In other words, applying the method 
to current data will not provide sufficiently tight constraints on 
£2 m , Q.a to be competitive with current cosmological methods. 
Ultimately therefore, to address systematic errors and to turn the 
proposed cosmic ruler into a useful cosmological tool, a large 
sample of homogenous systems is needed, e.g. high-redshift 
early-type galaxies with <x, At, and the locations of the images 
measured with high precision. From the sample, one may elimi- 
nate problematic systems, such as systems with a lot of external 
shear, or non-simple lenses. Moreover, each lens would have to 
be modeled in detail, i.e. accounting for systematic uncertain- 
ties by measuring them directly and including the effects in the 
analysis. 

Future facilities will allow such an experiment. LSST , 
SKA, JDEM, Euclid, and OMEGA (iDobke et alj 1200 



cal galaxies appear to be fairly universal dBak & Statlerll2000l 
lAlam & Rvdenll2002l) . The effects of extrinsic mass are impor- 



Ivezicetal.l l2008t ICarilli & Rawlingsl 12004 iMarshall et al. 
20051: iMoustakas et alj 12001 " will both find large numbers of 
new lenses and have the potential to accurately constrain the 
time delays. And with the James Webb Space Telescope (JWST) 
or future 30^40-m class telescopes (e.g. TMT - Thirty Meter 
Telescope and E-ELT - European Extremely Large Telescope), 
velocity dispersions can be obtained to the precision required 
and velocity anisotropy ca n be effectively cons trained through 
integral field spectroscopy dBarnabe et al.l 120091) . This must be 
coupled with high-resolution imaging which can constrain both 
the density structures of the lenses and significant mass struc- 
tures affecting the lensing. Extrinsic mass and density structures 
can also be constrained indirectly by modeling the lensing con- 
figuration, including flux ratios and positions of the images. 
The Dark Cosmology Centre is funded by the DNRF. 
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